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Demonstration of the stability or instability of 
multibreathers at low coupling 



Abstract 

Whereas there exists a mathematical proof for one-site breathers stability, and an 
unpublished one for two-sites breathers, the methods for determining the stability 
properties of multibreathers rely in numerical computation of the Floquet mul- 
tipliers or in the weak nonlinearity approximation leading to discrete non-linear 
Schrodinger equations. Here we present a set of multibreather stability theorems 
(MST) that provides with a simple method to determine multibreathers stability 
in Klein-Gordon systems. These theorems are based in the application of degen- 
erate perturbation theory to Aubry's band theory. We illustrate them with several 
examples. 
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1 Introduction 



Discrete breathers are time-periodic, localized oscillations in discrete systems 
due to a combination of nonlinearity and discreteness. They have become a 
well understood phenomenon since the publication of the proof of existence in 
Ref. [1]. This proof is the origin of exact and powerful numerical methods to 
calculate them and to determine their stability [2]. A deeper insight has been 
achieved since the introduction of Aubry's band theory [3]. The last reference, 
together with Ref. [4] can be considered as reviews, although a new one is 
badly needed after the huge development of the subject in the last years. 

The stability of one-site breathers is proofed in Refs. [2,3] under rather general 
conditions, the stability of the possible two-sites breathers in analyzed in an 
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unpublished theorem in Ref. [5], page 69. Other methods use the hypothesis of 
small amplitude oscillations or the rotating wave approximation [6,7] leading 
to the Nonlinear Schrodinger equation (DNLS), which is the actual equation 
analyzed, and, there are, of course, the efficient but slow numerical methods 
mentioned above. 

Here we propose a method based in the properties of the bands in Aubry's 
bands theory and perturbation theory, to obtain the stability properties of 
any multibreather at low coupling. In many cases it involves only some modi- 
fication of the coupling matrix and the knowledge of the hardness/softness of 
the on-site potential. We will briefly summarize the band theory in the next 
section, while we introduce the notation and some basic concepts, on which 
our theory is based. However, we refer to Ref. [3,1,5] for detailed explanations, 
as it will be long and repetitive to expose it in detail. 

In section 3 we develop the method for symmetric on-site potentials or in- 
phase multibreathers, which is synthesized in a theorem, commented in sec- 
tion 4. Its scope is enlarged to non-symmetric potentials in section 3, and 
to generalized Klein-Gordon systems in section 6. The method is applied to 
several interesting examples in sections 7 and 8. In Appendix A we calculate 
the value of a magnitude 7 used in our theory, and in Appendix B we relate 
the curvature of the bands with the characteristics of the on site potential. 



2 Band theory and notation 

2.1 The Newton operator 

We consider Klein-Gordon systems with linear coupling described by dynam- 
ical equations of the form: 

N 

u n + V'(u n ) + e C nm u n = n — l,...,N (1) 

m=l 

where the variables u n are functions of time t, V{u n ) is an homogeneous on- 
site potential, V its derivative, ii n denotes derivation with respect to time, 
N is the number of oscillators, C is a coupling constant matrix, which can 
describe nearest neighbor or long-range interaction, and includes the boundary 
conditions, and e is the coupling parameter. We use a notation similar to 
quantum mechanics, i.e., \u) = [ui(t), . . . , ujv(i)]^ meaning the transpose 
matrix). Defining V{u) = [V(ui), . . . , V(un)]^ and analogously its derivatives, 
equation (1) can be written as: 

+ \V'(u)) + eC\u) = 0. (2) 
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We suppose that the functions u n (t) are time-periodic (with period T and 
frequency u^,), time-reversible solutions and, therefore, they can be written as 
cosine Fourier series with real coefficients (some of these assumptions will be 
relaxed later). The (linear) stability properties of a given solution \u) depend 
on the properties of the characteristic equation for the Newton operator Af £ 
given by 

K(u) |0 = ||) + V"{u) ■ |0 + e C |0 = E |0 , (3) 
where • product is the list product, i.e., f(u) ■ |0 is the column matrix with 
elements f(u n (t)) £ n (t). If E — 0, this equation describes the evolution of small 
perturbations |0 of \u). 

2.2 The Floquet matrix 

Any solutions of Eq. (3) can be determined by the column matrix of the initial 
conditions for positions and momenta Q(0) = [£i(0), . . . £/v(0), 7Ti(0), . . . ttn(0)]\ 
with 7i n (t) = £ n (t), and Q(t) describes its evolution in the space of coordinates 
and momenta. A base of solutions is given by the 2 N functions with initial 
conditions ft v (0), v = 1, . . .,2 N, with ftj'(O) = 8 vh We will identify often a 
given solution £(£) of the Newton equation with the corresponding matrix of 
initial conditions Q(0). As the Newton operator depends on the T-periodic so- 
lution u(t), it is also periodic, and the evolution of the solutions of Eq. (3) can 
be studied by means of the Floquet operator, which maps £(0 into £(i + T). 
In a finite system this is equivalent to the Floquet matrix Te being given by: 



'{UT)} 




'{Uo)} 


K(T)}_ 




.K(o)}_ 



(4) 



These matrices can be easily constructed numerically, by integrating Eq. (3) 
2N times from t — to t — T, with initial conditions f^(0). Then, the v 
column of Te is given by Q U (T). 

2.3 Stability and bifurcations 

The eigenvalues {Xi}f£[ of Fq, called the Floquet multipliers, determine the 
linear stability of the solution u(t). If there is any eigenvalue with |A;| > 1, 
the corresponding eigenfunction 0(0 grows with time and u(t) is unstable; if 
A; = 1, 0(0 is T-periodic; if A; = —1, 0(0 is 2 T-periodic. However, F (and 
any Fe) is a symplectic matrix, as it is derived from the symplectic system 
Eq. (3), which implies that if A; is a non-zero eigenvalue, so it is 1/A. Fe is also 
real, which in turn implies that Af and 1/A* are also multipliers. Therefore, 
the eigenvalues of F E come in groups of four, if they are not real, or in pairs, if 



3 



they are real or if | A; | = 1. We can write the Floquet multipliers as {exp(i#/)}, 
with 0i, in general, complex numbers. The complex number \9i are called the 
Floquet exponents, and 0\ the Floquet arguments. 

The only possibility for the stability of u(t) is that all the eigenvalues of F 
have moduli 1, i.e., they are at the unit circle. Therefore, the condition for 
linear stability of u(t) can be described as all the Floquet arguments for Fq 
being real. Moreover, if u(t) is stable and, then, every |Aj| = 1, the Floquet 
multipliers come in complex conjugate pairs with arguments (±0j) or as double 
1 or — 1 (0 = or = ±7r). If a parameter like the coupling e is changed the 
multipliers of F change continuously. Therefore, a bifurcation to an unstable 
solution can only take place in three different forms: a) Harmonic instability: 
two complex eigenvalues moving along the unit circle collide at A = 1 (0 — 0); 
b) Subharmonic instability: two complex eigenvalues moving along the unit 
circle collide at A = — 1 (0 — ±7r); c) Oscillatory or Hopf instability: two pairs 
of complex eigenvalues moving along the unit circle collide at ±9 ^ and 
abandon the unit circle as a quadruplet (A, 1/A, A*, 1/A*). 

A Floquet multiplier of F is always known. Calculating the derivative of 
Eq. (2) with respect to time we obtain that M ii = 0. u is also T-periodic, 
therefore, it is an eigenfunction of Fq with multiplier 1, and as they come in 
pairs, there is always a double 1 multiplier, ii is called the phase mode because 
its meaning is that if ii(t) is a solution of Eq. (2), u(t + dt) is also a solution. 
While the solution u(t) exists, this double eigenvalue is always there. Due to 
the possible forms of the bifurcations, if all the multipliers are isolated except 
the double 1, the system is structurally stable, i.e., there is a neighborhood in 
the space of the parameters where u(t) is stable. 



2.4 Aubry's band theory 

However, it turns out that much information can be obtained by studying also 
the Floquet arguments for E 7^ 0, with is known as Aubry's band theory [3]. 
The set of points (9,E), with a real Floquet argument of F E have a band 
structure. The fact that the Floquet multipliers come in pairs of complex 
conjugate pairs brings about that if (0, E) belongs to a band, (—0, E) does it 
too, i.e., the bands are symmetric with respect to and dE/ d9(0) = 0. As u 
has eigenvalue 1 (0 = 0), there is always a band tangent to the axis E = at 
9 = 0. There are at most 2 N points for a given value of E and, therefore, there 
are at most 2 N bands crossing any horizontal axes in the space of coordinates 
(9,E). The condition for linear stability of u(t) is equivalent to the existence 
of 2 N bands crossing the axis E — (including tangent points with their 
multiplicity). If a parameter like the coupling e changes, the bands evolve 
continuously, and they can loss crossing points with E = 0, bringing about 



4 



the instability. 



We know more about the eigenfunctions of the Newton operator. The fact 
that it is periodic, or, in other words, that it commutes with the operator V of 
translation in time a period, Vf(t) = f(t + T), implies that its eigenfunctions 
can be chosen simultaneously as eigenfunctions of V, which by Bloch theo- 
rem, are given by £(9,t) = x(9,t) exp(i Ot/T), x(9,t) being a column matrix 
of T-periodic functions. It is straightforward to check that exp(i#) is the cor- 
responding Floquet multiplier, and 6 its Floquet argument (Fe is simply the 
representation of V in the base 

2.5 Bands at the anticontinuous limit 

The key concept on the demonstration of breather existence [1] , single breather 
stability [3] and the present paper is the anticontinuous limit, i.e., the system 
with all the oscillators uncoupled, e = 0, in Eq. (1-2). At the anticontinuous 
limit, Eq. (1) reduces to N identical equations: 

u n + V'(u n ) = . (5) 

Supposing that we consider time-reversible solutions around a single minimum 
of V, there are only three different solutions: a) oscillators at rest u n = 0; 
b) excited oscillators with identical u n (t), hereafter denoted u°(t); c) excited 
oscillators with a phase difference of n with the previous ones, given by u n (t) = 
u°(t + T/2). Each site index is given a code a n , which takes elements in 
{0,1,-1}, where o~ n = represents an oscillator at rest, u n = 0; o~ n = 1, 
an oscillator with solution u°(t); and o~ n = —1, the solution u°(t + T/2). The 
matrix of codes a = [a±, . . . ,ctn]^ represents the state of the system at the 
anticontinuous limit. 

We suppose that there are p oscillators at rest and N — p excited oscillators. 
Equation (3) at e = becomes 

^o(«) |0 = 10 + V"(u) ■ |0 = E |0 , (6) 
or, equivalently, N — p identical equations: 

'L + V'\u n )i n = E^n, (7) 

with only a periodic eigenfunction for E = 0, the phase mode u n {t) for the iso- 
lated, excited oscillators. The other one is the growth mode, (see Appendix B) 
which is not bounded and will not be used in this work. They give rise to N—p 
bands tangent to the axis E — 0, shown in Fig. 1. They have positive cur- 
vature at (6, E) = (0, 0) if the on-site potential is soft ( dH/ du>b < 0) and 
negative if it is hard (dH/ du b > 0), as it is demonstrated in Appendix B. 
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The possibility of dH/ cM, = is excluded by the conditions of the breather 
existence theorem [1], which means that the on-site potential is truly nonlin- 
ear for the isolated oscillators with the breather frequency u^. If, as we are 
supposing here, the on-site potentials V(u n ) are identical, the solutions u n (t) 
are also identical except for a change of phase and the bands are superposed, 
if not, we might have different bands but with the same general shape. 

The remaining p equations corresponding to the oscillators at rest are of the 
form: 

L + M 2 £n= E^n, (8) 

with ujq = \JV"(0). They have only a T-periodic solution for E = , the null 
solution, due to the non-resonance condition for breather existence nub ^ u)$ 
(p G N) [1], i.e., none of the harmonics of the breather resonates with the oscil- 
lators at rest. They provide p identical bands which are easily calculated solv- 
ing the equation above. Its solutions are £ = exp(±iy / cjQ — E t), with Floquet 

multipliers exp(±i\J ljq — E T) and Floquet arguments 9 = ±\Jljq — E T. 
That is, the bands are given by E = — cu^(9/T) 2 , where 9 can be re- 
duced to the first Brillouin zone [— 7r, 7r] by the addition of 27rp, p e Z. If the 
oscillators have different rest frequencies, the rest bands are not superposed 
but they have the same characteristics. 

The rest bands are also shown in Fig. 1. Note that this figure is only a sketch, 
for clarity, as very often some excited oscillators bands are very flat and dif- 
ficult to appreciate at the same scale. This sketch reproduces, however, the 
basics facts of the band structure. 

When the coupling is switched on, the degeneracy of the bands is generically 
raised. All the tangent bands at (9, E) = (0, 0) except one, than continues 
there, move upwards or downwards. If the on-site potential is soft, Fig. 1 
(left), and a band moves upwards, a double tangent point with the axis E — 
is lost, or, in other terms, a pair of Floquet arguments of Fq becomes complex 
and the solution u(t) is unstable. If the on site potential is hard, Fig. 1 (right), 
the same occurs when a band moves downwards. 



3 Multibreathers stability 

In this work use degenerate perturbation theory [8] to demonstrate the stabil- 
ity or instability of the breathers of any code. Degenerate perturbation theory 
establishes that if A/o is a linear operator with a degenerate eigenvalue E , 
with eigenvectors {|u n )}, which are ortonormal with respect to a scalar prod- 
uct, i.e., (v n \v m ) = 5 nm , and if e H is a perturbation of A/o, with e small; then, 
to first order in e, the eigenvalues of Ao + e N are E +e Aj, with Aj being the 
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eigenvalues of the perturbation matrix Q with elements Q nm = (v n \N\v m ). 
Note that perturbation theory as described in the reference cited, is time inde- 
pendent perturbation theory, but our time variable t is their spatial coordinate 
x. 

Let us consider again Eq. (3) with zero coupling, e = 0. As explained in the 
previous section, if there are N — p excited oscillators, there N—p zero eigen- 
values corresponding to the T-periodic eigenfunction u n (the phase modes of 
the isolated oscillators), or in other words, a N—p times degenerate eigenvalue 
Eq = if we restrict the domain of Af to periodic functions. What we need to 
know is the sign of this degenerate eigenvalues when the coupling is switched 
on. 

To apply degenerate theory to multibreather stability we need to identify the 
perturbation operator, a suitable scalar product and a ortonormal basis of the 
eigenspace with eigenvalue E = 0. The scalar product is defined as: 

N r T/2 

&I6> = E / -df. (9) 

n=l J - T l 2 

Let us suppose initially that all the excited oscillator are identical and vibrate 
in phase. We will denote by u°(t) these identical solutions. In this case, all the 
phase modes are also identical and will be denoted as ii . The N — p elements 
of the basis are 

\n) = i[0,...,0,M°,0,...,0]t, (10) 
H 

with the non-zero element at the position n, n being the index of the excited 
oscillators, and /i = \J J^ 2 / 2 (u ) 2 dt. It is straightforward to check that they 
are ortonormal. 

It is enough to know the eigenvalues corresponding to periodic and real so- 
lutions, i.e., with Floquet argument 6 = 0, as the intersections of the bands 
with the axis 9 = correspond to periodic solutions of (3). Therefore the 
N — p |n)'s form the basis of the degenerate eigenvalue E = needed to 
apply perturbation theory. 

To obtain the perturbation operator, we expand in Taylor series Eq. (3) at 
e = and obtain to first order in e : 



KH |0 = \0 + V"{u) ■ |0 + e (V"'{u) ■ u £ ■ |0 + C |0 ) = 

= (^ ) + eA i )|0> (11) 



with u £ = (§sf) and u is also the e = solution. 



s=0 
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Fig. 1. Sketch of the band structure for a soft on-site potential (left) and a hard 
on-site potential (right). The continuous lines correspond to excited oscillators and 
the dashed ones to oscillators at rest. 

Therefore the perturbation operator of Ao is 

M\i) =V"'(u)).u £ .\0 + C\0, (12) 

where u and u £ are calculated at e — 0. We do not know u £ , but deriving with 
respect to e, at e — 0, the dynamical equations (2) we obtain: 

\u e ) + V"(u) ■ \u £ ) + C\u) = or N \u £ ) = -C\u) . (13) 

The perturbed eigenvalues Ei are e \ , Aj being the eigenvalues of the pertur- 
bation matrix Q, with elements (n\J\f\m). The matrix C of elements {n\C\m) 
is simply the matrix C without the columns and rows corresponding to the 
oscillators at rest. The other terms are: 



(n\V"'(u) ■ u £ \m) = 




with u Hj£ = y-ffjj^j ■ Thus, only the diagonal elements (n\V"'(u) ■ u £ ■ \n) 
are non-zero. To calculate the last integral in (14) we will integrate by parts 
and use that the integral in a period of the derivative of a periodic function is 



8 



zero. Besides, the functions u nj£ are periodic as the coefficients of their Fourier 
series are given by the derivatives with respect to e of the Fourier coefficients 
of u n . In the deduction below, all the integral limits are —T/2 and T/2, and 
the terms between brackets from integration by parts will be zero. The last 
integral in Eq. (14) becomes: 



u°u n , £ V"(u°) 



V'(u u ) ii r , 



T/2 

-T/2 
T/2 

-T/2 



J V"(u°)u°u n , £ dt - I V"(u°)u°u n , £ dt = 
+ J V\u )u £<n , dt- j V"(u°)u°u nt£ dt = 

- / U°(Un,e + V"(u°) Un , £ )dt (15) 



The term between parentheses, is the n component of the lhs of Eq. (13), 
i.e., it becomes -E m CnmMm) where v? m = u°, if the oscillator m is excited, 
and zero otherwise, i.e., it is — I] m C nm -u° = — (J2m C n m) u°. Equation (15) 
becomes: 



(E^J / u°u°dt = 

m 

(J2C nm ) ([u°u ] T/ * - f{uydt) = -(£C nm V 2 . (16) 

That is, equation (14), leads to: 

(n\V'"(u)-u £ -\n) = -J2C nrn at e = 0. (17) 

m 

Therefore the diagonal elements of the perturbation matrix Q are 

Qnn ^ ' C nm ~\- C-an ^ ' C nm ^ ] Qnm ■ (^^) 

To summarize, the perturbation matrix Q is given by: 

Qnm C nm , 71 ^ Tfl , Qnn ^ ' Qnm 5 (^^) 

Vm^n 

C being the coupling matrix without the p rows and columns corresponding 
to oscillators at rest. 

This result can be very easily extended, to the case when there are oscillators 
out of phase, at least for an even potential V(u). In this case u°(t + T/2) = 
—u°(t) and ii°(t + T/2) = — u°{t). Let us define a new code a, which is equal 
to a, but with +1 instead of zeros, and perform the ansatz u n = a n u n . By 
substitution in the dynamical equations (1), and taking into account that 
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V'(u n ) = V'(a n u n ) = a n V'(u n ), because V'(u) is odd, we obtain: 

a n u n + a n V'{u n ) + C n . (20) 

m 

Multiplying by a n , and using that its square is +1 we get: 

u n + V'(u n ) + VnC nm CT m U m = 0. (21) 
m 

As the elements of C with an index corresponding to an oscillator at rest will 
not appear in the calculations, we can define a new coupling matrix, with 
elements (? n C nm a m) i.e., Diag(cr) CDiag(o-). The eigenvalues Ei will have the 
right sign, but the eigenfunctions £ of the Newton operator corresponding to 
(21) have to be transformed to recover the original ones by |£) = Diag(a) |£) 

The result can be summarized in the following way: 

Theorem 1 (Symmetric MST) Consider a Klein- Gordon system, with ho- 
mogeneous on-site potential V , linear coupling matrix C and coupling param- 
eter e, and a multibreather given by a vector code a at zero coupling, with p 
zero elements, and N — p non-zero elements. The potential V must be even 
if there is any —1 in a. We define the following matrices: S = Diag(a); C 
the N — p, squared matrix equal to S C S but without the rows and columns 
corresponding to the zero codes in a; Q, the perturbation matrix, with non- 
diagonal elements Q nm = C nm , and diagonal elements Q nn = — J2\/ m ^n C nm . 
Let {Aj}^ P be the eigenvalues of Q and suppose that there is only one zero 
among them. Then: 

a) ifV(u n ) is hard and there is any negative value in {e Aj} the multibreather 
at low coupling will be unstable, and stable otherwise. 

b) ifV{u n ) is soft and there is any positive value in {e\i\ the multibreather 
at low coupling will be unstable, and stable otherwise. 

There is always a zero eigenvalue is explained below. If there are more than 
one the stability is undefined within first order perturbation theory. 



4 Comments 

4-1 Global phase mode 

The perturbation matrix Q has always a zero eigenvalue corresponding to 
the global phase mode. It can be useful to make it apparent as it reduces the 
order of the secular equation by an unity. Considering the matrix Z — Q — XI , 
det(Z) = is the secular equation. Its roots are the same if we change Z to 
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Z' obtained by summing to the first column the rest of them, i.e., Z' nl = 
- Evm^n C n , m - A + Evm^n C'n.m = -A. Now, subtracting the first row from 
the rest of them, the first column has only a non zero element, which is equal 
to A. Therefore, A = is an eigenvalue of Q and the secular equation is reduced 
by an unity. 

The corresponding eigenvector x satisfies Q x = 0, i.e., — (Evm^n C nm ) x n + 
X^Vm^n C nm x m = 0, Vra, with solution x n = 1, Vn. Therefore, x = S x = a, 
which corresponds to the global phase mode. 



4-2 Extended and reduced forms 

Note that the code a is the matrix of components of the multibreather at zero 
coupling in a N-dimensional basis, { \n)}, with the only non zero element u° 
at the position n, including the indexes of all the oscillators. For numerical 
calculation and interpretation of the eigenvalues and eigenvectors of Q it can 
be more convenient to define C = S C S, i.e, conserving the rows and columns 
corresponding to the rest oscillators but all their elements set to zero. It will 
be denoted the extended form while the previous one is the reduced form. The 
price for using the extended form is to obtain p extra zero eigenvalues, which 
causes no harm, except if there are more than p + 1, as it can make more 
difficult the interpretation of the eigenvectors. 



4-3 Diagonal elements 

The diagonal elements of C play no role, because they do not appear in the 
coupling matrix. Therefore a coupling of the spring-type as e ( 2 u n — u n -i — 
u n+ i) gives identical results as the dipole-type — e (w„_i + u n+ i). 



44 2D and 3D lattices 

Note that the theorem applies equally to 2D and 3D lattices, provided that 
we can number the sites with an index n, or consider n as a multi-index, and 
construct the coupling matrix C. 
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5 Non-symmetric potentials 



The extension to non-symmetric potentials V(u n ) is straightforward, but it 
involves some numerical or analytical calculation. That is the reason for not 
including it in theorem 1, and we present it here as another theorem instead. 
If the potential is non-symmetric, there are two time-symmetric, T-periodic 
solutions to the equations for the isolated oscillators at zero coupling, Eq. (5): 
u°(t), with code +1, and u°(t + T/2), which is no longer —u°(t), with code —1. 
The N —p elements of the basis, \n) will be as in Eq. (10) but with u°(t + T/2) 
instead of u°(t) if a n = —1. Let us define the symmetry coefficient, 7, in the 
following way: 

J T { 2 /2 u°(t)u°(t + T/2) dt 
_ tU^Al K — — (22) 

G /2 u°(t)u°(t) dt 

If {z k }l™ are the Fourier coefficients of a A; m -truncated Fourier series of u°(t), 
i.e., u°(t) = z Q + Efe=i2-2 fe cos( ku b t), then: 

7 = T k m ,0 2 • (23) 

Note that 7 depends on the breather frequency through the {z k }, and it is 
positive and smaller than 1, being its value 7 = 1 for the symmetric potential 
case. Following the deduction of theorem 1 it is straightforward to obtain: 

Theorem 2 (Non-symmetric MST) The generalization of theorem 1 for 
non-symmetric potentials consists of constructing the perturbation matrix Q 
in the following way: 

a) If n or m are indexes corresponding to oscillators at rest Q nm = 0. 

b) If n and m, n 7^ m, are indexes corresponding to excited oscillators in 
phase, i.e., o~ n o~ m 1 ; Qnm C nm . 

c) If n and m, n 7^ m, are indexes corresponding to excited oscillators out of 
phase , i.e., a n a m = - 1, Q nm = -^C nm . 

d) The diagonal elements are Q nn = —J2\/ m ^nQnm 

e) The rest of the procedure is exactly the same as described in theorem 1, 
getting rid of the rest columns and rows (reduced form) or not (extended form), 
calculating the eigenvalues {\i}i~ p of Q and obtaining the stability properties 
of the multibreather from the signs of the {e\i} and the hardness/softness of 
the on-site potential. 

The coefficient 7 depends on the on-site potential and the frequency u h . We 
demonstrate in the Appendix A that 7 = oo h for the Morse potential. Figure 2 
shows some numerically calculated values for different potentials. They tend 
to +1 when uj^ — > u = 1 because at this limit u°(t) = A sin(cc>bt) and 
u°{t + T/2) = —u°(t) as in the symmetric on-site potential case. Generally 
speaking, if the frequencies are not very far away from the linear frequency 



12 



1 

0.9 
0.8 
0.7 
0.6 
f-0.5 
0.4 
0.3 
0.2 
0.1 








// 






/ / 






/ / 


- 




/ / 






/ / 






// 






1/ 






1/ 


- 




/l 






/ 1 


- 




/ 1 
1 

1 

1 




/ / 
/ / 
/ / 

/ / 






/ / 






/ / 






/ / 






/ / 






/ 






/ 






/ 






/ 









0.2 



0.4 0.6 
CO, 



0.8 




Fig. 2. Some values of the symmetry coefficient 7. Left, soft on-site potentials: 
dashed line: V(u n ) = ^u^ — \u^\ continuous line: V(u n ) = ^(exp(— u n ) — l) 2 . 
Right, hard on-site potential: V{u n ) = + \v? n + \u\. Rest frequency ojq = 1. 

u>o, 7 is usually close enough to the unity and the sign of the eigenvalues do 
not change with respect to the symmetric on-site potential case 7 = 1. 



6 Generalization 



Theorems 1 and 2 are equally applicable if the system has inhomogeneities at 
the coupling matrix but not if they are at the on-site potential. Here we intro- 
duce the generalization to a broad class of Klein-Gordon systems, including 
different masses or inertia momenta, inhomogeneous on-site potentials with 
several minima, nonlinear coupling, solutions with different phases (i.e, not 
restricted to time-reversible solutions). Let us consider a Klein-Gordon sys- 
tem with Hamiltonian 

H = E (\ m n < + V n (uS) + e W{u) , (24) 

and dynamical equations 

x dW(u) . , 

m n u n + V;(u n ) + e—^- = 0, (25) 
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which at zero coupling become: 



m n u n + V^{u n ) = Vra. (26) 

Suppose we can label the solutions of these equations, with given frequency uj\, 
and period T, with a code composed of an index for the site n, another index 
for the different minima and solutions (biperiodic, etc, different potential wells) 
a n , which can be different for each site. If we considering a periodic solution for 
the n-th equation , it can be written as a 2 7r periodic function g n u) h a n (^ht+a). 
Let us suppose that g n uj h a n is chosen such that a = corresponds to the time 
symmetric solution with g nu)ha „(0) = and g nUb<Tn (0) > g n u h a n {^)- Therefore, 
the solutions of equations (26) are determined by the code {n, a n ,a n }. If we 
consider a given solution at the anticontinuous limit we can simply write it 
as u = [ui, . . . , u° N }\ but now, the are, in general, different, because of the 
particular multibreather solution u, whose stability we are interested in. 

The characteristic equations for the Newton operator Af £ are given by: 

. d 2 W(u) 

(A4 £)n = mnU + KTKKn + £ £ 7j U = E £ n (27) 

which to first order in e calculated at e = lead to 

(A/q On + e (A? On = (E + e A) £„ , (28) 

where: 



(A/o On = m n Cn + K'K) C„ = So £n , (29) 

d 2 W(v°) 

(A^On = «„, £ + £ o o Cm = A£ n . (30) 



with u Hi£ = (^~) • The first set of equations is the characteristic equa- 
tion for the Newton operator at zero coupling, and the second one for the 
perturbation operator Af. The only periodic solution of each equation (29) 
with Eq = for a given solution u of (26) , is Thus, we can construct a 
ortonormal basis of ker(iVo), restricted to periodic solutions (and, therefore, 
that can be chosen real) of given period T, as \n) = — [0, . . . , 0, u° n1 0, . . . , 0]*, 

with the only nonzero element at the n position, and fi n = \j j^^i^n) 2 
Deriving Eq. (25) with respect to e at e = we get: 

(AT -) n = m n u n , e + = - £ (3D 
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Following exactly the same procedure as in section 3 we get the perturbation 
matrix Q, in extended form, with non diagonal elements for indexes corre- 
sponding to excited oscillators: 



Q 



nm 




u 



W{u°)\m) 



J-T/2 




U 



"n 



d 2 W{u°) 
du n du m 



u m dt , n ^ m (32) 



and diagonal elements 



Qnn ^ ' Q 



nm ■ 



(33) 



If there are p oscillators at rest and N — p excited oscillators, the extended 
matrix Q has p zero columns and rows corresponding to the indexes of the 
oscillators at rest. The reduced form of Q is the N — p squared matrix, equal 
to Q but stripped of those zero rows and columns. As a result we can establish 
the theorem: 

Theorem 3 (Generalized MST) Given a generalized Klein-Gordon sys- 
tem (24), a specific multibreather solution at zero coupling {u^}, determined 
by a suitable set of codes, u the corresponding solution at low coupling, {X}^J[ P 
the eigenvalues of the reduced matrix Q, with only one zero, and being the wells 
corresponding to the N — p specific w° , of the same type, hard or soft, then: 
The solution u is stable if: 

a) The on-site potentials are soft and there is not any positive value in {e Aj}^ ; 

b) The on-site potentials are hard and there is not any negative value in {e Aj}^ 

This theorem includes both theorems 1 and 2. Certainly, it is of not so straight- 
forward applicability, although for certain cases there are only a few matrix 
elements Q nm to calculate. As an example, if there are just two types of sites, 
the on-site potentials are symmetric, the coupling is linear and we consider 
only time-reversible multibreathers, there are only two different matrix ele- 
ments to calculate. 



7 Applications 

7.1 Two-site breathers 

A trivial example is the theorem about the stability of the two-site breather 
found in Ref. [5] page 69. Here we generalize it to non-symmetric on-site 
potentials. The coupling matrix corresponding to periodic boundary condi- 
tions and of the standard (spring-type) attractive form, has elements C nm = 
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2<5 nm — £|n-m|,l an( i CVl — Cin — 



— 1, that is: 



2 - 



1 ... 



-1 



-1 



2-1 ... 







C = 



(34) 



-1 



... 



-1 



2 



If the two oscillator are in phase, the perturbation matrix is 



1 -1 



Q 



-l 



l 



with eigenvalues (Ai, A2) = (0, 2), the first one corresponds to the phase mode, 
the second to the antisymmetric eigenmode. Thus, one of the two bands tan- 
gent to the axis E — at zero coupling will remain there, and the other will 
move upwards. If the on-site potential is soft, i.e., with positive curvature, a 
intersection/tangent point will be lost, and the two site breather is unstable. 
For the hard on-site potential, the curvature of the bands is negative, and 
when the band moved upwards, the intersection points will be kept, and the 
system is be stable. If e < the signs of the eigenvalues e are reversed and 
so are the conclusions. 

For the out-of-phase, two-site breather, with code a = [1,-1], the perturba- 
tion matrix is: 



with eigenvalues (Ai, A2) = (0, —27). The conclusions above are reversed, being 
the instability mode, when appropriate, the symmetric one. Note that the 
boundary conditions in C have no effect. 

Theorem 4 (Two-site breathers) Consider time-reversible, two-site breathers 
in an homogeneous Klein-Gordon system with attractive linear coupling, then: 

a) The two-site breathers with codes ± [1, 1] and soft on-site potential are un- 
stable and with ± [— 1, 1] are stable. 

b) The two-site breathers with code ± [1, 1] and hard on-site potential are sta- 
ble and with ± [—1, 1] are unstable. 

c) If the coupling is repulsive, the conclusions are reversed. 



Q = 



-7 7 



7 -7 
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7.2 Three-site breathers 



Let us renumber the sites 1, 2, 3 and suppose that the code of the site in the 
middle is a 2 ,— +1 as a reference, i.e., a = [ai,l,a 3 ]. We define a, = 1 if 
<7j = 1 and a, L = —7 if a { = —1, for i = 1,2. Using C from Eq. (34) the 
perturbation matrix is: 



Q 



a\ —a\ 
-Ct\ Ct\ + «3 —0:3 
-a 3 a 3 



(35) 



Its nonzero eigenvalues are \± = ct\ + a 3 ± \J {a\ + o^) 2 — 3 «i 0:3. In order 
to have both eigenvalues with the same sign, we need that at\ 0:3 > 0, i.e., either 
«3 = ai = 1, or a 3 = ai = —7. That is, if <j\ = o" 3 = 1, \± = 1,3 > 0, if 
<j\ = S3 = —7, A-t = —7, —37 < 0, and if a± a-i — — 1, then A± have different 
signs. Therefore: 



Theorem 5 (Three— site breathers) Consider time-reversible, three-site breathers 
in an homogeneous K-G system with nearest-neighbor, linear coupling, attrac- 
tive for e > 0. Then: 

a) The three-site breathers with codes a = ±[—1,1,-1] are stable if e > 
(e < 0)and the on-site potential is soft (hard). 

b) The three-site breathers with code a = ±[1, 1, 1] are stable if e > (e < 0) 
and the on-site potential is hard (soft). 

c) All other three-site breathers are unstable for any sign of e and type of 
on-site potential. 

We can compare the predictions of the values given by the theory with the 
numerically calculated values of the intersection points of the bands with the 
axis 6 = 0. As an example, the numerically calculated eigenvalues and the 
predictions given by theorems 1 and 2 for the three-site breather with code 
[— 1, 1, —1] are plotted in Figure 3. The very good accuracy is evident. 



7.3 A model with long-range interaction 



An example of system with long range interaction is the twist model used in 
Ref. [9]. The dynamical equations are: 

T w, x , "^ /2 cos[^(n-m)] 
u n + V (u n ) +e 2^ j- _ u m = (36) 



m=n-N/2 



\n — m\ 
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Fig. 3. Comparison between the eigenvalues E of the Newton operator obtained 
numerically (points), predicted by the symmetric MST, theorem 1, (dashed lines) 
and by the non-symmetric MST, theorem 2, (continuous line), which is the right 
one to be used. Morse potential, code [—1, 1, —1], frequency lo^ = 7 = 0.8. 

V being the Morse potential, and 6 tw the angle between two dipole moments 
corresponding to neighboring base pairs in a simplified model of DNA. The 
coupling matrix elements are: 

_ cos[9 tw (n - m)] 
\n — m\ 6 

Some of the results, also checked numerically and with good accuracy with 
respect to the eigenvalues E, given by theorem 1, in spite of V being non- 
symmetric are: 



Code 


@tw — 


Qtw = 7T 


11 


Stable 


Unstable 


1 -1 


Unstable 


Stable 


101 


Stable 


Stable 


111 


Stable 


Unstable 


1-1 1 


Unstable 


Stable 


11-1 


Unstable 


Stable 



As an example of the effect of the symmetry coefficient 7 in theorem 2, we can 
consider the code [1, —1, 1] and plot the Newton eigenvalues E as a function 
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Fig. 4. Eigenvalues A of the perturbation matrix Q for a three-site breather with 
code [—1,1,-1] on the twist model. Continuous line: 7 = = 1; dashed line: 
7 = = 0.8; and dotted line: 7 = ui^ = 0.1. Only values of 7 very far from 1 lead 
to a different stability prediction. 

of 6 tw . The perturbation matrix is: 



7a — b —7a 6 
—7a 27a —7a 
6 —7a 70 — 6 



(38) 



with a = cos(# tU) ) and b = cos(2# to )/8. Figure 4 shows the dependence of the 
eigenvalues A on the twist angle 9 tw for three values of the symmetry coefficient 
7. It can be seen that only values of 7 very far from 1 lead to a change on the 
stability prediction. 



8 Multibreathers, phonobreathers and dark breathers. Parity in- 
stabilities. 



For larger, time-reversible multibreathers in homogeneous Klein-Gordon sys- 
tems with nearest-neighbor, linear coupling, it is trivial to obtain the stability 
properties of the two different solutions with wave number q — 0, all the os- 
cillators in phase, or q = ir, neighboring oscillators out of phase. Consider 
any number N' of contiguous oscillators in phase in a system coupled by the 
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matrix C in Eq. (34). The N' x N' perturbation matrix is given by: 



Q = 



(39) 



1-10... 
-1 2-1 . . . 

0... -1 2-1 
... 0-1 1 



The characteristic equation for Q x = A x corresponds to the following equa- 
tions: 



X\ — X2 = A X\ , 

2x n — x n+ i — x n ^i = Xx n , Vn 1,2, 

x jv' — xn'-i = \ xn> ■ (40) 

Substitution of the trial solution x n = cos(gn — <f> q ), leads to 4> q = q/2, q = 
mn/N', with m — 0, . . . , N' - 1, and A m = 4 sin 2 (m 7r/(2 iV')), i.e., positive 
eigenvalues without a degenerate zero. In case all the oscillators of the group 
are out of phase, Q — > — and A m — > — 7A m . Therefore: 

Theorem 6 If e > 0, the time-reversible, in-phase, multibreathers are un- 
stable (stable) with soft (hard) on-site potential. For any out-of-phase multi- 
breathers, and e < the conclusions are reversed. 

Although, we have no mathematical proof, according to numerical calculations 
of the eigenvalues of the perturbation matrices corresponding to groups with 
different codes, the numbers of negative and positive eigenvalues are equal to 
the numbers of —1 and +1 in {a n c^+i}^! 1 - Therefore being unstable for any 
on-site potential. 

Note that if there are several non-contiguous multibreathers, the perturbation 
matrix if composed of independent blocks of perturbation submatrices of the 
same type, each one with a zero eigenvalue by construction as seen in subsec- 
tion 4.1. The stability of such a system cannot be determined by the MST. 
The only possible conclusion is that if there are eigenvalues of different sign 
this system would be unstable. 

Phonobreathers are multibreathers without oscillators at rest. For time-reversible 
phonobreathers in a system either with free-ends, or fixed-ends boundary 
conditions, the perturbation matrix Q is the same as in Eq. (39). If the 
boundary conditions are periodic, the only change is that Qn — Qn,n = 2, 
Qni — Qi,n = — 1, and the eigenvalues are the same. That is, the mode q = 
with attractive coupling and hard potential is stable. Changing the hardness, 
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the type of coupling, or the mode, changes each time the stability. In this 
way, we reproduce the results in Refs. [7,6] as a consequence of modulational 
instability and using the DNLS approximation, respectively, and in theorem 9 
in Ref. [3], based in the properties of the action. 

Parity instabilities are a interesting effect for phonobreathers with periodic 
boundary conditions. It is a consequence of the parity of the total number of 
sites N. If, for example, V is soft, the coupling is attractive, and N is even, 
all the eigenvalues of the perturbation matrix are negative, i.e., the system 
is stable. But if N is odd, there is an isolated positive eigenvalue, bringing 
about the instability. There is a clear physical origin of this phenomenon, as 
the boundaries are equivalent to an inhomogeneity. According to our theory 
it is also obvious: we have, in fact two oscillators at the ends coupled and in 
phase, which gives rise to a positive eigenvalue. 

Dark breathers are mult ibreat hers with only one or a few oscillators at rest. 
Their stability depends on the characteristics of the contiguous groups of ex- 
cited oscillators, as a consequence of theorem 6. We have to be careful about 
the number of particles, for example the it mode with code around the dark 
site [. . . , 1, —1, 0, 1, —1, . . .], needs an odd number of sites to avoid having 
two in-phase oscillators at the ends, and a parity instability. The one with 
[. . . , —1, 1, 0, 1, —1, . . . ] needs an even number of sites. If the system has free 
of fixed periodic conditions, the perturbation matrix is decoupled in two and 
the stability is undefined as commented above. 

We have compared the predictions of the instabilities and the values of the 
intersection points of the bands with the axis 9 = 0, in many cases, as for 
example in Ref. [10], with very good results. The question of up to which 
values of the coupling parameter e the theory is valid, is however unsolved. 
Generically speaking values up to e = 0.1 (compared to the rest frequency 
coo — 1), are safe, but there are some exceptions. In the same reference, for 
the dark breather hard on-site potential and attractive coupling, a mode not 
predicted by the MST gives rise to an instability as soon as e = 0.022. 

As the width of the phonon band can be calculated, and the width of the 
background is given by the eigenvalues of the perturbation matrices, it is also 
possible to predict the values of the e, for which oscillatory and subharmonic 
instabilities occur, but we do not extend further here. 



9 Summary and conclusions 

We have developed using degenerate perturbation theory a method for ob- 
taining the stability properties of mult ibreat hers of any code at low coupling. 
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The method is synthesized in three different versions of a theorem, referred to 
as the multibreather stability theorem or MST for short. The two simplest ver- 
sions correspond to time-reversible multibreathers, with linear coupling and 
homogeneous on-site potentials. If the on-site potential is symmetric or we 
consider only in-phase multibreathers, it involves only some simple modifica- 
tions of the coupling matrix. If none of those conditions are fulfilled then it 
involves the analytical or numerical calculations of a magnitude 7, called the 
symmetry coefficient, as its value is 1 for symmetric on-site potentials. For 
soft potentials and values of the frequency until about one half of the rest 
frequency, a good approximation is the multibreather frequency, result which 
we demonstrate is exact for the Morse potential. For hard potentials its value 
its close to unity. The generalized version of the MST is of not so straightfor- 
ward applicability, although its complexity depends on the characteristics of 
the system and multibreather considered. 

We give some examples of application of the method and compare it with nu- 
merical results. The systems considered are two and three-site breathers, with 
nearest neighbor or long-range interactions, multibreathers, phonobreathers 
and dark breathers. A parity instability can appear in finite systems depending 
on the parity of the number of oscillators. All these examples are interesting in 
themselves, but also illustrate the use of the MST, and show that it provides 
a powerful method of easy applicability to determine multibreathers stability. 

Some other applications under study are multibreathers in 2D and 3D systems, 
prediction of the oscillatory and subharmonic instabilities and the peculiarities 
of some disordered systems 
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A Calculation of the symmetry parameter 7 for a Morse potential 

In this appendix, we calculate the value of the symmetry parameter 7 in 
theorem 2 for the particular case of a Morse potential. The choice of this kind 
of potential relies in the fact that the expressions of the orbits are easy to 
manage. We will proof that, for this potential, 7 = u^. 

Let us suppose an isolated oscillator submitted to a Morse potential. The 
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energy of the system is: 

£ = \x 2 + \{e W {-x)-l?- (A.l) 

We look for solutions with period T = 2ix/uo h that are time-reversible, i.e, with 
x(0) = 0, or, in other words, at t — the system is at a turning point. There 
are two of them, obtained from the equation above: X\ = — log(l + \/2E) < 
and x 2 = — log(l — \f2£) > 0. If we consider the solution with x(0) = xi, t(x) 
is given by: 

t(x) = / , A.2 

^2(5 - (exp(-rr) - l) 2 /2) 

which leads to: 

1 ■ r 1 (, 25-1 m\ 
t(x) = ; arcsm -= (1 + ) . (A.3) 



VI - 2£ \2 L V2£ V exp(-x) 

By substitution of t — T/2 = 7r/cJb and x = x 2 we obtain c^b = \/l — 25. 
The inversion of Eq. (A.2) leads to: 



1 — a/1 — WuCOSOA.^ 
X(0 = log * ^ • ( A - 4 ) 



To calculate 7, we need the Fourier coefficients of x(t). If the solution x(t) is 
expressed as: 

00 

x(t) = z Q + 2 2 fe cos(fca;b , (A-5) 



fc=i 



the coefficients {z k }^ =0 are given by: 

1 r T / 2 

Zk = — / cos(A;a;b £) dt . (A.6) 

T J-T/2 

The calculation of these integrals is straightforward and leads to: 

z = log , 2 ; z k = --\ 1. 2 ) ■ ( A - 7 ) 

b 



2 ^b n\Jl_ UJ 2 



The parameter 7 is given by: 

7 = - E f ( ~iiy - (a.8) 

By substitution of the Fourier coefficients in Eq. (A.7), we obtain 

7 = - E ^ lH f with r = l^. (A.9) 
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As |r| < 1, 7 can be easily calculated using the expressions for the sum of an 
infinite geometric progression: 



leading to 7 = c^b, as we wanted to proof. 



B Curvature of the bands 



In this section we demonstrate a key point used in our theory: the fact that the 
band corresponding to an isolated excited oscillators has negative curvature 
if the on site potential is hard and positive if it is soft. Unfortunately, the 
demonstration is quite long if we want it to be self-consistent. A similar one 
can be found in Ref. [11] and different ones in Ref. [5,12]. 

Let us consider the dynamical equation of an isolated oscillator with Hamil- 
tonian H = 1/2 p 2 + V(u), with u a single-value, real function of time and 
p — it: 

u + V'(u)=0 (B.l) 
The characteristic equation for the Newton operator Af corresponding to a 
given periodic solution of Eq. (B.l)is 

Aft = Z + V"(u)Z = EZ, (B.2) 

where £ is a C 2 function of t. For each eigenvalue E there are only 2 independent 
solutions, determined, for example, by the values of £(0) and £(0). We can 
obtain them for E = 0, by deriving Eq. (B.l) with respect to t and with 
respect to Uy,: 

u + V"(u)u = <=> Nil = ; u w + V"(u)u w = <=> Nu w = . (B.3) 

As we have a single oscillator we can suppose that u(t) is time symmetric with 
a suitable origin of time, i.e., u is a time-antisymmetric function of t with the 
same period T = 2 7r/c<Jb as u(t): u(t + T) = u(t), u(0) = u{T) = 0. The 
properties of u w = du/dui,, can be also obtained. We can write u as a cosine 
Fourier series u(t) = J2k ) =o c k cos ( UJ bt), with the coefficients depending on 
uj h . Therefore, du/du h = J2kLo(dck/duj h ) cos(cJ b t) + J2kLo ~ ^ * s m (^b t), or: 

u w = -^- = 7 (i) + —u(t) (B.4) 

uCc>b 

7(t) being a time symmetric, periodic function of time with the same period 
T as u. This functions constitute a base of the eigenspace corresponding to 
E = 0. 
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On the other hand, Af (B.2) is a periodic operator because it depends on a 
periodic function u. Its eigenfunctions are given by the Bloch theorem: 

Z(e,t)= X (W,t)exp(&t/T), (B.5) 

with x a T-periodic function of time and 9 = 9(E). If we consider only values 
of E for which 9 is real, the set of points (9, E(9) constitutes a band. If £ is 
a solution of Af £ = E£, the fact that Af is a real operator implies that £* is 
another solution, i.e., (—6,E) also belongs to the band and this is symmetric 
with respect to 9 and, therefore, dE/ d9(0) = 0. As Af is time-symmetric 
£(6, —t) is another, and has to be proportional to £* as there can be only two 
independent solutions. 

The Floquet matrix Fe is defined as: 





= F E 




_£(T)_ 







(B.6) 



The eigenvalues {A/} of Fe are the Floquet multipliers, if we write them as 
\i = expi#z, with 9 real or complex, i6i are the Floquet exponents, and 9i, the 
Floquet arguments. An eigenfunction with 9i real is bounded, if 9\ = 1 the 
eigenfunction has period T. The condition for linear stability of the solution 
u is that all the Floquet arguments of F are real. 

It is easy to check some properties. The Bloch functions Eq. (B.5) are eigen- 
functions of Fe with Floquet exponent 9, therefore Fe is diagonalizable (over 
C) with Floquet multipliers exp(±i#) and Floquet arguments ±9. However, 
at E — 0, there is only one eigenfunction of F , u, the other independent 
function of the subspace Af £ — 0, u w , it is not an F -eigenfunction. That is 
F has a degenerate Floquet multiplier 1, or Floquet argument and can be 
transformed into a Jordan block. Another consequence of the latter is that 
(9, E) = (0, 0) belongs to the band. We need to calculate the curvature of the 
band at this point. 

Let us come back to the Bloch functions £(9,t) in Eq. (B.5) and £*(#,£) 
corresponding to a pair of symmetric points (±#, E) of the band, and sup- 
pose that 9 — > (and then E — > 0). At 9 = the two functions collide in 
£(0, t) = x((M) — X*(0,t), therefore £(0,t) is a real, T-periodic function. But 
we already now that at E — 0, there is a T-periodic function solution of 
Aft; = 0, u. Therefore, £(0, t) = u (Af is lineal, so we can adjust the norm of 
£ in order that the latter equation is fulfilled). 

We can reobtain the missing function that spans ker(A/"). If we derive Eq. (B.2) 
with respect to 9, we obtain Af £e — E e ^ + E^ e and at E = (E e (0) = 0) we 
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obtain: jVf fl (0,t) = 0. Deriving f(0,t) in Eq. (B.5) we get: 
& = Xe(0, t) exp(i£/T) + i^ X (9, t) exp(i£/T) = X e exp(i£/T) + i^ ■ (B.7) 

At 9 — 0, £ = u and has to be proportional at u w , which we know is the 
other independent function. Therefore 

. * t ui h { T t \ luI ( 2tt t \ 

T T Vicjb w b / 27r y uj h J 

Comparing with Eq. (B.4) we obtain 

^(0,t) = i^u w . (B.8) 

The following step is to express the curvature Egg of the band E{9) at (9, E) = 
(0, 0) in terms of u and its derivatives. Let us derive Eq. (B.2) twice with 
respect to 9, we obtain: 

Uiee = Eeei + 2E e ie + Eiee^Miee = Eeeii at 9 = . (B.9) 

Multiplying to the right by u and integrating a period we get 

E ee [ T u 2 dt= [ T iiNieg&t. (B.10) 
Jt=o Jt=o 

For any pair of functions f(t) and g(t) 

f fHg = [ T (fg + fV"g)dt=[fg]J =0 -[fg]J =0 + f (gf + gV" g) dt 
Jt=o Jt=o Jt=o 

= [fg}J=o-lfg}J=o+ [ T g^ffat. 

Jt=o 

Applying the latter equation: 

[ T uMfa dt = [uiee]J =0 - [u&e]I =0 + / iteeMu) dt . (B.ll) 
Jt=o Jt=o 

Using that J\fii — 0, u(0) = u{T) = and u(0) = u(T) the previous equation 
leads to 

( T uM&edt=-u(T)(&e(0,T)-&e(0,0)) . (B.12) 
Jt=o 

To calculate the last difference, let us derive Eq. (B.7) with respect to 9. We 
obtain 
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Zee = Xeo exp(i6t/T) + Xo exp(i0t/T) + & 

= xoo exp(i0f /T) + i| (6 - i| + 4 6 

t t 2 
= X^exp(i^/T) + 2i-£ e + — £ 

* t 2 

=>• fee = Xoo + 2i- 6 + j^u at = 0. 

X and therefore xeo are T-periodic functions, therefore £,qo{0, 0) = Xoo(®, 0) 
X^(0, T) and &*(0, T) - ^(0, 0) = 2i &(0, T) . Then 



I 



T r)H ; 2 ; 2 H 

«jve w dt = -ii(2i&(o,r)) = ^n 2i (i/^m) = ---3- • 

Substituting in Eq. (B.10), and naming the action I = J^qU 2 dt we finally 
obtain at 9 = 0: 

M°) = -4^- CB.I3) 
Til du h 

I is always positive, therefore, if the potential is hard (dH / duj^ > 0), the 
curvature Egg is negative and the band is tangent to the axis E — from 
below, if it is soft ( dH/ du\, < 0) the curvature is positive and the band is 
tangent form above, as we wanted to demonstrate. 

Note that this demonstration is also valid for the whole coupled system with 
the obvious changes: u — > [ui, . . . , wat]^; V(u) — > V(u) + eW(u), W(u) being 
the coupling potential; V'(u) — > [dV/du\, dV/du n] V"(u) the matrix 
with elements d 2 V/du n du m and analogously. We also need the additional con- 
straint that the solution u(t) is chosen time-symmetric and that the pair 
of eigenvalues with 9 = is isolated. This means that the band tangent at 
(9, E) = (0, 0) continues having the same curvature properties depending on 
the hardness/softness of the potential 



References 



[1] RS MacKay and S Aubry. Proof of existence of breathers for time-reversible or 
Hamiltonian networks of weakly coupled oscillators. Nonlinearity, 7:1623-1643, 
1994. 

[2] JL Marin and S Aubry. Breathers in nonlinear lattices: Numerical methods 
based on the anti-integrability concept. In L Vazquez, L Streit, and VM Perez- 
Garcia, editors, Nonlinear Klein-Gordon and Schrddinger Systems: Theory 
and Applications, pages 317-323. World Scientific, Singapore and Philadelphia, 
1995. 



27 



[3] S Aubry. Breathers in nonlinear lattices: Existence, linear stability and 
quantization. Physica D, 103:201-250, 1997. 

[4] S Flach and CR Willis. Discrete breathers. Physics Reports, 295:181-264, 1998. 

[5] JL Marin. Intrinsic Localised Modes in Nonlinear Lattices. PhD thesis, 
University of Zaragoza, Department of Condensed Matter, June 1997. 

[6] AM Morgante, M Johansson, G Kopidakis, and S Aubry. Standing waves 
instabilities in a chain of nonlinear coupled oscillators. Phys. D, 162:53, 2002. 

[7] I. Daumont, T. Dauxois, and M. Peyrard. Modulational instability: fisrt step 
towards energy localization in nonlinear lattices. Nonlinearity, 10:617-630, 
1997. 

[8] D Griffiths. Introduction to Quantum Mechanics. Prentice Hall, New Jersey, 
USA, 1995. 

[9] B Sanchez-Rey, JFR Archilla, F Palmero, and FR Romero. Breathers in a 
system with helicity and dipole interaction. Phys. Rev. E, 66:017601-017604, 
2002. 

[10] A Alvarez, JFR Archilla, J Cuevas, and FR Romero. Dark breathers in Klein- 
Gordon lattices, band analysis of their stability properties. New Journal of 
Physics, 4:72.1-72.19, 2002. 

[11] T Cretegny. Dynamique collective et localisation de I'energie dans les reseaux 
non-lineaires. PhD thesis, Ecole Normale Superieure de Lyon, 1998. 

[12] G Kopidakis and S Aubry. Intraband discrete breathers in disordered nonlinear 
systems I: Derealization. Physica D, 130:155-186, 1999. 



28 



